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I. INTRODUCTION 



Until high energy experiments achieve the neccesary energy and sensitivity to discover 
new particles in direct production, any experimental progress beyond the Standard Model 
(SM) should go through the observation of loop effects in precision measurements. These 
new effects can be parametrized using effective Lagrangians where coefficients of higher 
dimensional operators carry all the information on new physics. In this paper I ask how 
accurately one can measure these couplings in proton-proton collisions at presently available 
energies. 

The heavy mass of the top quark allows us to speculate that new physics could show 
up more easily in the couplings of the top quark then of any other SM particle. The 
most obvious place do such precision measurements is the Fermilab Tevatron, where at 
Run 2 an estimated 20,000 top pairs will be produced. In this paper I am looking for small 
new operators, generated in many theories with strong electroweak physics, large extra 
dimensions, etc. The succes of the SM leads us to expect that these contributions are small, 
and this in turn allows us to keep only the interference with the leading SM process, top 
pair production through qq — > g — > it. The emerging restriction singles out eight different 
operators and these can be disentengled once we make use of the information residing in the 
spin state of the produced top quarks. This spin information is in turn translated into the 
angular distribution of the decay products of the two t, t quarks. 

The effect of new physics operators in top pair production has been repeatedly studied 
in the literature Jj]-|5[]. Most of these investigations look at top transverse momentum and 
invariant it distributions in addition to effects in the total cross section. In the following 
I argue that this approach looses most of the information present in the events because 
it sums over top polarization states. The top quark decays into t —>■ Wb so quickly that 
its spin state is not washed out by strong interactions. In Refs. the use of particular 
asymmetries helps to incorporate at least some of this information. We will see however, 
that the full use of angular distributions significantly increases the accuracy. 

The first attempt in this direction was taken in 0, where the effect of the top chromo- 
magnetic moment on these distributions was analyzed. We now extend this analysis for all 
operators that can contribute, up to dimension 6. In fact, there is little reason, other than 
theoretical prejudice based on adherence to particular sorts of models, to expect that any 
of the operators in the set should be distinctly larger than others. For example, the dimen- 
sion 5 chromomagnetic moment K is produced at O(-p), similarly to dimension 6 operators 
(A is the scale of the new physics.) In this paper I purposefully avoid any reference to the 
expected relative size of the new couplings. 

In what follows I classify the observable contributions in two cases: (i) when the elec- 
troweak (EW) symmetry breaking sector contains only a light Higgs doublet, and (ii) when 
all new fields are heavy. In the former case we can write an effective Lagrangian where the 
fields are in a linear representation of the EW group, while in the latter case we must use a 
nonlinear representation. The approach is thus quite general and should only break down if 
there are new resonances in the below-f TeV region that interact significantly with the top 
quark (such as top pions, for example.) We write down the new operators in both cases and 
conclude that this experiment cannot tell apart the two scenarios. 

The paper is organized as follows. In Sec. |I| I first discuss the relationship of the two 
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representations, and write down a basis of those new operators that contain at least one 
top field. Then I find the corrections to the git vertex, to top decay in the tWb vertex and 
4-quark vertices: two CP violating and six CP conserving operators. I proceed in Sec. |TTT| by 
calculating the contribution of each of these new interactions to the it spin density matrices 
and to top decay, and compute their contribution to the fully differential cross section in 
the 6-particle phase space. These contributions can be used for comparing experimental 
distributions to disentangle each contribution. 

It is important to note that one does not need to reconstruct the probability of events in 
the many-dimensional phase space from the measurement. That would be impossible with 
the event numbers involved (on the order of several thousand). One can instead estimate 
the new couplings, using for example a Maximum Likelihood method. We perform such 
an analysis in Sec. [TV], using Monte- Carlo- generated data according to the SM distribution. 
This analysis tells us the attainable statistical accuracy of the estimation. 

I sum the results in the Conclusion and discuss the discovery potential of TeVatron Run 2 
at Fermilab. I find that new physics at a scale order TeV can be generically detected and 
that all-hadronic decay modes of the top quarks may play an important role, the complexity 
of these events notwithstanding. 

II. LINEAR AND NONLINEAR REPRESENTATIONS: PARAMETRIZATION 
OF NEW PHYSICS IN INTERFERENCE 

We are interested in effective Lagrangians controlling top pair production and subsequent 
decays with or without a single SM Higgs field, where all other physics resides in the higher 
dimensional operators. Without including any model dependent "prejudice", we should 
proceed by writing down all possible operators at each dimension level and look for their 
experimental consequences. 

When a SM Higgs boson is present in the Lagrangian, the underlying SU{2) x U(l) 
symmetry imposes restrictions on the possible form of the new operators. We are free to 
choose the fields in a linear representation of the gauge group and implement spontaneous 
electroweak breaking by introducing a Higgs vacuum expectation value. The introduction 
of a Higss v.e.v however induces lower dimensional new operators of a non-gauge invariant 
form. 

It may well happen that the SM Higgs is very heavy or is not present at all. In the 
absence of light resonances an effective Lagrangian still make sense, but in that case its 
field content does not allow us to put all the fields in closed linear representations. Because 
we are supposing that all additional fields are heavy, we are allowed to put the SM fields, 
including the three eaten-up Goldstone bosons, in a nonlinear representation of the gauge 
group and require gauge invariance of the effective Lagrangian. 

It has been shown H, however, that the additional restrictions by imposing this nonlin- 
early realized symmetry are exactly compensated by the additional freedom in the choice of 
the gauge. Note that the presence of unphysical degrees of freedom in a generic gauge (such 
as longitudinal and timelike components of the gauge bosons) lets us write down additional 
operators containing these fields. In the unitary gauge all these nonphysical fields fall away 
and the net result is that all restrictions on the new operators from gauge invariance dis- 
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FIG. 1. The only types of new top physics processes that contribute to qq — > tt in interference. 

appear. The gauge non-invariant operators allowed this way will be similar to those in the 
linear representation that contain no Higgs fields. 

In this section I write down, in both cases, all the new physics operators that can effect 
the qq — > tt process, up to dimension 6. There is a huge number of such operators. However, 
one may impose the following restrictions on those that may give appreciable contributions: 

• The use of an effective Lagrangian presupposes that the corrections are small, therefore 
we are justified to keep only the interference between new physics operators and SM 



• We keep new physics only in vertices that involve at least one top quark field. Other 
new physics, including that in light quark interactions, can presumably be seen easier 
in other processes. 

• We drop those new operators that pick up a weak coupling from an unmodified SM 
vertex (such as corrections to the Ztt vertex, contributing to qq —>■ Z — > tt). 

• We drop all loop corrections. 

With the above criteria some combinatorics shows that the only contributions come from 
(a) corrections to the git vertex, (b) chirality preserving 4-quark operators of the color octet 



interfere with the SM process, and (c) corrections to the tWb vertex in top decay (see Fig. [[]). 

The linear set of 2-quark (plus boson fields) operators have been classified extensively 
in the literature 0||,[nj. Any complete prediction of the contributions makes sense with 
reference to a complete set of operators, because any partial list will tacitly imply a particular 
choice of operators that have zero contribution. For this reason I provide in Appendix [A] 
a list of the operator basis I am using. This list completely agrees with that of Ref. |9| 
when all CP-odd operators are dropped. In Appendix [B] I also provide the corresponding 
operators in the nonlinear set. 



processes. 




where the light quarks q, q are must be of the same flavor to 
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The contributions we find to the gtt vertex are only: 



ys 4m t 2 ys 4m t 



'* ta^—*G a ul/ t, 



(i) 



where the nonlinear set contributes k = /t 33 and k 



/t 33 , while the linear set contributes 



k = —4m t -^Rebll and k = Am t ^Imfr 3 ]] (for the detailed meaning of the r.h.s. see the 
Appendix.) Both k and k are required to be real by Hermiticity. The absence of the 0(p 2 ) 
terms in the expansion of the gtt form factors seems rather strange. These operators have 
been turned into 4-quark operators by the use of the equations of motion. 

Next we look at corrections to top decay. We find the following three operators that 
seem to contribute to the tWb vertex: 
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^= b 3 , 3 (linear) 



-%^W»J R a^b L + h.c. 
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c 3 
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(2) 



In the following we will see however, that only three real linear combinations of these op- 
erators contribute to the differential cross section; the operator with c 3 does not contribute 
at all. It is not unexpected to see fewer operators contribute to top decay because of the 
kinematical restrictions. Only the CP-even combination 



A*i = 2 

and the CP-odd combination 

A*2 = 2 



ml - 
ml 



ml - m^ 



ml 



Re (2c 2 + ci + c 3 ) 



Im (2c 2 + ci + c 3 ) 



(3) 



(4) 



affect the differential cross sections, while the CP-even 

^ = -2 m "-f w lm(c 1 + c 3 ) 

ml 



(5) 



only changes the total rate by an overall factor of (1 + (Iq). Because this last factor would 
show up only in the total top decay rate, it is irrelevant for our investigation and I drop it in 
the following. Note that in a similar manner I ignored the contribution of a ~ W 1 1^ 7^ b^ 
operator which would only modify V t b- 
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Finally, we turn to the effect of 4-quark operators. These are the same in both set and 
are severely restricted by the above criteria: 



9s f ( ^ \ (- ,A a \ ( A a \ /_ „ ,A a 




4m? x | C yy x ^ y gj ^ y tj + c VA x ^ y gj ^Vy*) (6) 
+cav x ^7„7 5 y^ f^y*) + caa x ^7 

where the coefficients cyy, c va, c^y, caa, must all be real (q = u, d, s, ... is the light quark). 
The selection of these operators has been pressed on us by the required presence of interfer- 
ence with the SM process. However, our operators can be embedded into the set of custodial 
invariant 4-fermion operators discussed in [[J. Therefore, we find no direct restriction on 
our operators from custodial symmetry. I separated off a factor of g 2 s for easy comparison 
(this factor does not need to reflect any strong interaction physics.) 

This completes the list of all contributing operators. Their size remains remarkably un- 
restricted by previous experiments and by theoretical arguments. None of them contributes 
linearly at tree level to the electroweak precision parameters. Consistency of the effective 

theory requires that their size cannot significantly exceed (j^j , but the new physics scale 
A is unknown. All we know for sure is that they must be less than unity, and this fact 
justifies throwing away all contributions that are quadratic in the new couplings. 

It is now a straightforward but tedious calculation to compute the contribution of each 
of these operators to the differential cross section, a calculation to which we now turn. 



III. NEW PHYSICS CONTRIBUTIONS TO THE DIFFERENTIAL CROSS 

SECTIONS 

In this section I derive the most important result of this paper: the contribution of each 
of the operators in ([T],[2||J) to the fully differential cross section comprising ti pair production 
and decay. 

The SM process is shown in Fig. [| By D I designate the down-type quark in the decay. 
In (semi)leptonic decays it should be replaced by a charged lepton. In the following I will 
keep, for the economy of the notation, for the lepton momentum; I well proceed similarly 
with D, U and U. 

A Breit-Wigner formula is applicable to the amplitudes when the t and t, as well as 
the W ± are on shell: 

This equation allows us to separate the physics into top production, described by the 
density matrix Pss^-g and the decay represented by T®, and Tip: 

7f x p s - SjS ,s> = Y,\ M gg->tt\ 2 , (8) 
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FIG. 2. The full cascade of ti production and decay in the Standard Model. The decay products 
of the W + are an up-type U = u,c quark (or a neutrino) and a down-type antiquark D = d,s (or 
a positively charged lepton, e + or 



averaged for the incoming spin and color, summed for outgoing color (but the t,t spins are 
kept fixed). 

The decay of the t,t quarks is described by 
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In the decay contribution I separate out factors that do not receive contributions from 
new physics operators, 
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(10) 
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The six particle phase space becomes manageable when we describe the momenta of the 
three t decay products in the t frame and the three t decay products in the t frame. Denoting 
the three-direction of the W(W) by w(w) and the three-direction of the D(D) by n(n)in 
these frames, we find for the differential cross section 



[m 2 - m 2 w ) 2 



9 x 2 2 %9 m%m 2 w T 2 T 2 w (m 2 + m 2 w Y 
v(zmf) w Jf v„\ „ dQ(n) 



4 



9 



x 2 



m 2 + 



x \V u75 \ 2 \V VD \ 2 



xdn ( t> x dfi(n) x 



f3 w n ■ w) 



(Pb • Pu) x dft(n) x 
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(11) 
(Pb-Pu) x z t 



where dQ.[ ZMF ^ is the solid angle measure of the top direction in the zero momentum frame 

2 2 

(i.e. the center of mass frame of the tt system), (3w = ^fc^r is the speed of the W in the 
t frame, and the new physics corrections are all factorized into 

Zt = £ ^f- x r3 x r ®. (12) 

ss,s's 

The spin states corresponding to the indices s, s require careful definition: because we 
are also using the off-diagonal elements of the density matrix, the relative phase of the spin 
states must be taken into account. We achieve this using the following conventions. The 
spin vectors s,s lay in the scattering plane and their spatial components are back-to back 
to each other in the ZMF frame. A vector £ M is chosen normal to the scattering plane, and 
the for a spin basis I use the so-called "off-diagonal" basis of Ref. |TTJ in which the SM 



contribution to the density matrix is simple. Our conventions are best understood in the 
ZMF where the various vectors are 

p£ = ^(1|1|0|0)" pjf = ^(1| - 1|0|0)^ 

$ = cos 9|/3 sin e|0)^ V - t = ^(1| -/3cos6| -/3sin9|0)' 1 

^ = A iH 2n (lsinO|0|l|0)^ r = -r -J. 2n (/?sin9|0| - 1|0)" 

yl-/3 i sm / © y 1-/3 2 sin 2 © 

-((3 cos 9|1 - f3 2 sin 2 9|/3 2 sin 9 cos 9|0)^ 



sin 2 6 

7 

y/l-0 2 sir 
^ = (0|0|0|lf 



= . 7 . , (/3cos9| - (1 -/3 2 sin 2 9) | - /3 2 sin 9 cos 9 1 Of (13) 

•y 1— fP sin 



Here, (3 is the speed of the top quark in the ZMF frame, 7 = 1/yT — (3 2 and 9 is the top 
scattering angle. 

With these definitions the phases of the spin states are fixed by requiring 



M T M I = +\{l-Pt + m t) 7 5 7-(£ + iv) M l M t = +\il-Pt + m t) 7 5 7-(£ ~ 
^t 17 ! = -|(7-P? - mt) 7 5 7-(£ - "7) W = -|(7-P? - ™t) 7 5 7-(£ + *>7) 



(14) 



where the j (J,) states now have spin projection +1/2 (—1/2) on the corresponding spin 
vector. 

In the t(t) frame the purely spatial vectors r), s, £ (77, s, £) form right handed orthonormal 
bases, so that they are used as the respective X^y^Z £IX6S to help define the directions of the 
momenta of the decay products. Although the differential cross section in the six-particle 
phase space is independent of the choice of the t,t spin basis, the choice of the variables 
does depend on these Cartesian frames. 

Given these conventions we can now calculate the p and r matrices in the SM, 

Pn!n = Pmt = 2 - f sin 2 9, = ffi u = f sin 2 9, (15) 

~TT Tl 
IT II 



and all other elements vanish flTT|. The decay contributions are, in the notation 
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Jt) 

r SM 



St) 

T SM 



1 + s ■ n n-(£ + ir}) 

n-(£ — irf) 1 — s ■ n 

1 — s • n n- (£ + iff) 

n- (£ — iff) 1 + s ■ n 



1 + cos d sin de^ 

sin de~ lip 1 — cos d 

1 — cos $ sin de 1 ^ 

sin •&e~ %9 1 + cos 



(16) 



Substituting (0) and (|TBp into (0) we find 

^(SM) _ ^ _|_ cos ^ cos ^ )(i _ _ p 2 sin 2 6) + - (3 2 sin 2 9 sin^ sin^ cos(<y9 — Tp). 



(17) 



In these expressions I introduced the polar angles in the "off-diagonal" basis. These are 
defined asQ 



n • s = cos d n ■ £ = cos tp sin $ 

n • s = cos t? n • £ = cos </? sin $ 

w • s = cos x w ' £ = cos sin x 

w • s = cos x W ■ £ = cos sin \ 



n ■ r] = sin 99 sin 
n • 77 = sin </? sin $ 
w ■ r) = sin sin % 
w ■ ry = sin sin \ 



(18) 



Any corrections to the gti vertex can only show up in p while new physics in top de- 
cay will affect only the r's. At this point the astute reader may observe that while the 
decay contributions affect the distributions in all the observed angles, the corrections to top 
production cannot change the distribution in w, w. This clear distinction is, however, lost 
when experimental realities are taken into account. The jets produced by the U, D quarks 
are indistinguishable and the required symmetrization brings back the w-dependence. 

Now look at the contributions from each of the twelve new physics operators. I quote 
here only the result of a tedious but straightforward calculation. The contribution of the 
chromoelectric and chromomagnetic moments into Z t is 

Z K = k x {2 (1 + cos $ cos $) + /5 2 7 sin O cos B (sin $ cos $ sin tp + sin $ cos $ sin Tp) j 

+k x Ipj sinO \Jl — (3 2 sin 2 6 (sin cos $ cos p> + sin cos cos Tp) j , (19) 



and the four-quark operators contribute 

7 (2 — (3 2 sin 2 6) cos 6 (cos + cos §) — sin 6 (sin sin ip + sin sin Tp) 

Z c = c VA x [jj = 

y/l - p 2 sin 2 6 

+caa x /?7 {27 cos O (1 + cos $ cos $) — sin 6 (sin $ cos sin ip + sin cos $ sin j 
+c Ay x 2 7 2 y/l - (3 2 sin 2 (cos + cos ) (20) 
+cyy x 2 7 2 |(1 + cos?9cos^)(l - i/5 2 sin 2 6) + ^/3 2 sin G sin sin cos(</> — </?)>. 



1 There is a sign difference in the definition of the angies (p, p that I am using here and the one 
used in ||. The difference shows up in the form of the contribution due to n. 
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These contributions, as I explained above, do not depend on the angles x> 0; X; 0- The 
decay corrections, however, do: 



— iii x j sin^sinxcos((/? — 0) + sin $ sin % cos (y? — 0) 

+ (cos , # + cos i9)(cos x + cos x) x (1 ~~ s * n2 ®) 
— -/3 2 sin 2 9 sin -i? sin x cos — 0) + sin-$ sin x cos(y? — </ 
+//2 - \ sin $ cos $ sin x sin(<^ — 0) + cost? sin $ sin xsin(^ — 0) (1 — -j3 2 sin 2 9) 
— -/3 2 sin 2 9 sin d sin $ (cosx — cosx) sin(y2 — Tp) 

+ sin $ cost? sin x sin (</> — 0) + sin $ cos $ sin x sin 
;/i x >j (1 + cos^cos^Xl — -/3 2 sin 2 9) + -f3 2 sin 2 9 sin , (9sin'(9cos(v9 — Tp) 



The contribution of the CP conserving fii in the last line is proportional to the SM contri- 
bution and only rescales the CKM matrix elements by a factor (1 + [i{). Cne can drop this 
contribution. 

The form of the contributions in Eqns. ( P~D| , PU , 2~I ) constitutes the main result of the 
present paper. They can be used in an experiment to disentengle the contributions from 
the various new physics operators. Integrating over the various angles except 9 as required 
by Eqns. ( pTT]) immediately convinces us that all the information on many of the operators 
resides in top polarization. Namely, the integration is equivalent to replacing 



Z t Z t = M - -/T sin 2 9J [1 + ^ + 2Y x c vv \ + 2k + 2/V cos 9 x c AA - (22) 

The only couplings that affect the distribution of the events in the t,t momenta at all are 
fix, cyv and caa- Out of these, caa does not contribute to the total cross section. 



IV. SENSITIVITIES 



The parton-level calculations in Sec. |TXT| should be confronted with experimental realities. 
At first sight one would think that the large inaccuracies related to the extraction of parton 
momenta from jet observables preclude a detailed mapping of the phase space. However, in 
order to establish the presence of the new physics operators one does not need to reconstruct 
the full angular distribution. The complicated structure of Eqns. ( P^|PD| , |2l] ) does not allow 
us to select a simple observable that would indicate the presence of new physics. Instead, 
one should look for the particular types of angular correlations that are similar to those in 
the above equations. This is most easily done by a Maximum Likelihood estimate. Following 
this strategy the angular correlations are not spoiled by uncertainties in the jet variables 
even as large as ±50% or by very large backgrounds 0. 

The events fall into three categories, according to the decay mode of the t, t: 
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• Dilepton events. In this case two neutrinos are produced and full event reconstruction 
is extremely hard. The angles we are investigating are probably impossible to extract. 

• Semileptonic events. These provide an opportunity to accurately measure the momen- 
tum of the charged lepton (i.e. pjy or pd), but on the hadronic side the two light quark 
jets are practically indistinguishable. One must symmetrize the predicted differential 
cross sections and compare the resulting quantity to the experiment. Fortunately, 
we will see that this procedure does not entail much information loss. However, the 
smaller leptonic top decay branching ratio results in smaller statistics compared to the 
all-hadronic decays. 

• All hadronic events. In this case the advantages of high statistics are somewhat com- 
pensated by the hardships of relating six jets to six partons. Nevertheless, two SVX 
tags may identify the b, b quarks and requiring that the both tops and W's are on shell 
leaves a chance for an approximate event reconstruction. 

One is faced, however, with more complicated ambiguities in this case. There are two 
pairs of indistinguishable light quarks jets plus the entire t side is hard to distinguish 
from the t side. The former problem does not introduce large uncertainties in the 
estimation of new physics parameters. As for the latter, we observe that the effect of 
t t replacement amounts to the replacement of a <-> 7r — a of all the angles in the 
above formulas for Z t . Consequently, some of the operators (namely, caa and cav) 
are antisymmetric under this replacement and cancel out from the symmetrized cross 
sections. The rest however, including both operators in top decay, are symmetric and 
survive unchanged. 

In order to assess the statistical uncertainties involved in the determination of the new 
physics parameters I performed a Maximum Likelihood estimation of the eight real param- 
eters ^ = jtfj, f/,2, K, K, cyy, cyA, cavi c av- 1 used a Monte Carlo generator to produce events 
according to the SM distribution at y/s = 2 TeV and built a Maximum Likelihood estimator 
that finds the most probable values of the /Xj's (the true values are obviously zero.)Q The fall 
of the probability by a factor of e _1//2 was interpreted as a "1 — cr" interval (in the following 
I will call its size two times the "accuracy".) In each case the "accuracies" agreed well with 
the ones found from a second-order expansion of the log-likelihood function, so that the use 
of a covariance matrix is justified: 

ln£(/i) = ln£ ma:r - -//jCy^-. (23) 

As a check of the reliability of these estimates, I observed that in each direction in the 
parameter space when only one \ii is kept nonzero, the likelihood quickly decreases as \ii is 
moved outside the "1 — cr" interval. The true value /i, = has always been consistent with 



2 This analysis, an extension to many parameters of the one we used in Q, leads to the "best 
estimator" z ]j M -^Zt(m). After the publication of || I realized that this is the same estimator as 
the one described in | j!2| , found through a different argument unrelated to Maximum Likelihood. 
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the "accuracy." In addition, using different event numbers (all above M > 1200), I observed 
that the accuracy is indeed proportional to TV" -1 / 2 . 
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TABLE I. The statistical accuracy of the Maximum Likelihood estimate of each parameter ji^. 
The numbers mean 100 times the half-width of the "1-sigma" interval. The 2 nd , 3 rd and 4 th columns 
tell if the operator contributes to the total cross section, whether it survives t*-*t symmetrization 
or the integration over all decay angles except (called "No Pol."). The 2k, 20k refer to the 
number of events considered. "All Symm." denotes symmetrization for D *-* D, U <-> U and t*-*1. 
The column in boldface applies directly to semileptonic events at Tevatron Run 2. 



The accuracies of the results of this ML estimate are shown in Table [I] for 2000 and 20,000 
events. The former is relevant for the semileptonic events at the Fermilab Tevatron, while 
the latter is a hypothetical number shown in order to expose the gain in the accuracy when 
more events are observed. These numbers immediately confirm that the indistinguishability 
of the light quark jets does not drastically deteriorate the extraction, contrary to what one 
might have naively expected. 

The intricate pattern of correlations induced by the contributions to Z t tells us that 
the much of the information on the new operators resides in the spins of the t,t. When 
all the angles related to top decay are integrated out and only the variables (3 and 6 are 
kept, a reduced amount of information is still available. However, the only operators whose 
contributions do not vanish in this case are k, caa and cw We see in Table | that in the case 
of caa and eyy, much information resides in the @, 0-distributions and these two operators 
can be observed without measuring the angular correlations of the decay products. In fact 
Cw represents only a momentum-dependent, 0(p 2 ) correction to the strong coupling of the 
top, so its contribution is proportional to that of the SM process times r y 2 . This different 
dependence on the total energy allows that we see this operator in the unpolarized it cross 
section. In the case of k, the accuracy drastically deteriorates when we integrate out the 
polarization information, in accordance with the findings in Ref. [|1^ . 

The difference in the patterns of angular correlations introduced by each operator leads to 
the fact that, with one exception, the principal axes of the covariance matrix approximately 
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coincide with the chosen operator basis (i.e. Cy is approximately diagonal.) This means 
that the angular correlations provide almost independent measurements of each operator 
in our basis. In the case of k and cyy-, there is an (9(10%) mixing and the combinations 
corresponding to the principal axes are shown in the last two rows of Table |. 

The accuracies attainable in this way, supposing that systematic uncertainties will not 
play a drastic role, are generically ±0.01 to ±0.1. This can be translated to an observable 
new physics scale of A = 0(6 — 20)m t , with an uncertain numerical factor of order unityQ 
The 4-quark operators can be measured with a precision of ±0.01 in semileptonic events, 
as far as only statistical inaccuracies are taken into account. Two of these, cav an d caa 
are theoretically disfavored in many models as they involve an axial vector current of light 
quarks. The other two, cyv and cva may however remain unsuppressed as they are related 
to vector currents of the light quarks. Of course, the actual suppression of each operator is 
a model-dependent question and should not be addressed in this paper. 

V. CONCLUSIONS AND PROSPECTS OF DISCOVERY 

In the above we found that four CP-conserving and four CP- violating new physics opera- 
tors can contribute to top pair production in interference with qq — > it: the chromomagnetic 
and chromoelectric moments k,k, two decay corrections Ho, t^i and the four current-current 
color octet 4-quark operators cvv,cva,cav,caa (averaged over the light quark flavors.) In 
deriving our set of contributing operators I used the equations of motion, so that some 
of the physics in the git couplings may actually show up in the four-quark operators (for 
example, G^ ll ,i^'j^D u b is eliminated.) All these operators arise in both the linear and 
the nonlinear (chiral) set, so the two cases cannot be told apart in this process. I cal- 
culated the contribution of each of these operators to the fully differential cross section 
da(qq — > it — > 6 fermions), at tree level, in the approximation when all t,i, W + , W~ are on 
shell. This is the dominant process of top pair production on the Tevatron at ^/s = 2 TeV. 
At higher energies (for example at the LHC) our analysis is not relevant, because other 
Standard Model processes become important. 

Using a maximum likelihood method I estimated the statistical inaccuracy of the mea- 
surement of this parameter set. We find that only caa and cyy can be meaningfully measured 
without using the information encoded in the directions of the momenta of the six decay 
products. In addition, the contribution of all operators except cyy to the total cross section 
vanishes. It is well known that the total cross section is not measured very precisely at the 
Tevatron, due to inaccurate knowledge of the luminosity. This fact can affect therefore only 
the extraction of cyv but not that of the other parameters. 

The statistical inaccuracies in the case of 2000 semileptonic events (feasible at the Teva- 
tron), generically ± O(0.01 — 0.1), correspond to new physics on the TeV scale. The actual 



3 The normalization of the couplings m is so chosen that all mass dimensions are removed at the 
scale of ~ 2m t . Most it pairs are produced with close to threshold energies, due to the increase of 
the parton distribution functions at small x, so that rrif is the only energy scale in the problem. 
However, there remains some 0(1) arbitrariness in the normalization. 
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assessment of the meaning of this accuracy should be performed in the context of partic- 
ular models, but in any case these accuracies are an order of magnitude better than those 
previously found without the use of spin information. Although these results seem to be 
quite robust against large backgrounds or inaccuracies in event reconstruction, it remains to 
be seen how large is the effect of asymmetric systematic distortions of the distributions by 
experimental factors. This important question should be addressed by detailed simulations 
of the experiments, which the present author is not equipped to perform. 

As it was observed in ||, when the event number is less than ~ 1000, the estimation 
becomes useless due to very large inaccuracies exceeding the above quoted unitarity bounds. 
The small branching ratio of the leptonic top decay renders dilepton events useless for our 
purposes. The importance of the event number to improve the statistics leads us to consider 
the abundant all-hadronic decays modes. We have shown that the price to pay is the loss of 
the measurement of cav an d caa, but the increased event number can help to improve the 
accuracy of the rest of the parameters. 
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APPENDIX A: LINEAR REPRESENTATION 

In the linear representation there are no dimension 4 or 5 operators. In the following I 
provide a list of those operators with two quark fields, at dimension six, which in the unitary 
gauge generate vertices involving at least one top quark (terms proportional with light or b 
quark masses, and also all lepton fields, have been dropped): 



bf 




(Al) 


bf 


x (qliD^eiD^ + h.c. 


(A2) 


bf 




(A3) 


°4 


x{r L lJ^qi){^iD^)+h.c. 


(A4) 


h 


x (®iD^) + h.c. 


(A5) 


b aP 


x (KlA) (® T eiD^)+h.c. 


(A6) 


bf 


x (g£?4e$ T ) ($$) + h.c. 


(A7) 


g'bf 


xqla^B^u^ e$ T + /i.c. 


(A8) 


gK 


xg2^4 T -^V« v df + h.c. 


(A9) 


g s Wo 




(A10) 


ua/3 

gWi 


x?2^4 Y v ^® + h - c - 


(All) 
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Here, a, (3 are flavor indices, q% are the left handed quark doublets, u R , d R are the right- 
handed singlet quarks, B^ u and V^ u are the hypercharge and SU{2) field strengths, 
is the gluon field strength, and e is the antisymmetric 2x2 tensor in SU(2) space. The 
covariant derivative includes all electroweak and gluon fields. The coupling constants, 
bi, are (not necessarily Hermitian) complex matrices. In the above expressions the equations 
of motion have been used. This is consistent if one keeps the 4-quark operators . In the 
unitary gauge, replacing $ — > (0, allows us to write down the new physics operators in 
a physical basis where Goldstone fields are absent. 

There is a restriction on the size of 67 from the fact that it gives mass to the quarks. The 

absence of fine tuning requires that this additional mass, -|- y^J , does not significantly 
exceed the physical quark masses. This coupling, however, will not give a contribution to 

our process. The coupling 64 contains a piece — -|- u%^W^a L whose effect is to modify 
the corresponding CKM matrix element V a p. Again, the absence of fine tuning imposes a 
restriction on the size of 64. There are obvious restrictions of the flavor changing operators 
in the list. 



APPENDIX B: NONLINEAR REPRESENTATION 



In the nonlinear representation we have the following operators with two quark fields, 
again dropping those that contain no top quark: 



Dimension 4 

9 



2cosfl 



AgL x Z^ul + \f Ld x Z^dZ + \f Rll < Z /( -<;, -'■'„/,, 



9 
V2 



\w P R*W»U R ^d R + h.C. 



.a/3 
K ZR; 



(Bl) 



and \wl * s absorbed into the CKM matrix element V a @. Here, each Xz is Hermitian 
but Xwr can be any complex matrix. 

Dimension 5 

The magnetic moments are 



K a PpF 4- K a Pn —fl a I K a @ 9 V 
, n-v er UV ' ry n 9s n Lx,,,, "t- o „„„ a ^ U 



— U (7 



Arrii 



K a 0_9_W 

4m 4 



(B2) 



the electric moments are found from these with the replacement k — > k^X^ —>■ *X /1U . 
The weak derivative couplings are 



fa/3 ( / /a/3 

2cos# 2m 4 L ^ fl 1 V2 V 4m t 



.'a/3 



(B3) 
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Again, Hermitian and n' z , k' lw , k' rw need not be. Because 

Kw and Kw are not independent, k w *W^ v URa^d,L = ikwW^URa^di, we can require 
that k w and k w be both Hermitian. 

• dimension 6 

Dropping all operators that involve, in addition to two quark fields, at least two more 
electroweak bosons, we find for both H = L,R: 

a/3 

_2 Phz x ^ r i ^ U P (q Z + dvZ ) (B4) 
2 cos 9 ml 

a/3 

9 Phw x u a H Y i D 4 ( d„W v + d u W,) (B5) 



\/2 m 



t 

9 9s PHZ v 7 —a A ri xv 



2 COS 7?t| 



2 



^^x^^^r 7 ,4+^- (B7) 

9 9s PhZ 7 a a * r i!v „,/3 /oo\ 
X Z » U H^T G a 7i/«jj (B8) 



2 cos 6? 

^f^ x ^T^X+fc-c (B9) 

Here, Phz,p'hz,P~'hz are Hermitian but Phw, Phw Phw need not be - 

Note that in the nonlinear representation the covariant derivative involves only the 
gluon and electromagnetic fields but not the W and Z. All combinations with these fields 
should be included as there is no restriction from EW gauge invariance. 
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